Holstein model and Peierls instability in ID boson-fermion lattice gases 
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We study an ultracold bose-fermi mixture in a one dimensional optical lattice. When boson 
atoms are heavier then fermion atoms the system is described by an adiabatic Holstein model, 
exhibiting a Peierls instability for commensurate fermion filling factors. A Bosonic density wave 
with a wavenumber of twice the Fermi wavenumber will appear in the quasi one-dimensional system. 



The recent realization of Bose-Fermi mixtures of ul- 
tracold atoms P, S 3 is very promising for study- 
ing various strong correlation phenomena, with Bose 
fields replacing the lattice phonons of classic condensed- 
matter models. Virtual exchange of the bose-field exci- 
tations were shown to induce fermion-fermion attractive 
interactions |E IE S leading to Cooper-Hke pairing 
0, Hol llH and enhancing the transition to fermion su- 
perfluidity. Furthermore, completely novel phenomena 
were predicted, such as the formation of boson-fermion 
composite fermionic particles ll^ and their sub- 

sequent pairing into quartets |l2| . 

One particularly interesting direction is the possibil- 
ity of mutual fermion-boson trapping in optical lattices 
opening the way to the realization 
of various discrete lattice models. For electrons con- 
fined to a quasi-one-dimensional electric conductor the 
coupling to bosonic phonons leads to the Peierls insta- 
bility towards a charge density wave (CDW) with wave 
number equal to twice the Fermi wavelength 2kp. The 
corresponding instability for a Bose-Fermi mixture in 
an harmonic trap was recently predicted |l8l |. In this 
article we study the Peierls instability in a mixture of 
light fermionic and heavy bosonic atoms in a quasi-one- 
dimensional optical lattice, shown to be described by an 
adiabatic Holstein model. The resulting CDW depicted 
schematically in Fig. 1 , consists of both a fermionic den- 
sity wave and a spatial modulation in the bosonic density, 
with twice the Fermi wave length. Fermionic atoms and 
bosonic modulations will either position in alternate sites 
(Fig. la) or in the same sites (Fig. lb) depending on the 
sign of the fermion-boson interactions. 

We consider a mixture of Nc spin polarized fermionic 
atoms and Na bosonic atoms confined to a quasi one- 
dimensional (QID) optical lattice with AI sites (sec Fig. 

For a strong optical field, one can expand boson- and 
fermion field operators in terms of the one-mode-per-site 
Wannier basis set thus obtaining the lowest Bloch- 
band Hubbard model. 
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FIG. 1: Peierls instability in a lattice fermi-bose mixture: (a) 
repulsive fermion-boson interactions (b) attractive fermion- 
boson interactions. The shaded part depicts the bosonic mean 
field density whereas filled circles denote fermionic atoms. 



where the operators Cj and aj 



annihilate a spin polar- 



ized fermion and a boson respectively, in the j-th site. 
The density operators fij — c^Cj and n° = ajoj are the 
fermionic and bosonic densities respectively. The fermion 
and boson atomic masses and hopping amplitudes are 
mc,tc and ma-,ta respectively. The collisional terms gaa 
and gac correspond to on-site boson-boson interaction 
which will be assumed positive (repulsive) throughout 
the paper, and onsite fermion-boson interaction, respec- 
tively. The last term on the r.h.s. of Eq. ^ is the har- 
monic trap potential where luq, is the relevant oscillator 
frequency and £ is the lattice spacing. 



In the weak-interaction limit the bosonic field in Eq. 
may be treated in a mean-field approximation, re- 
placing the bosonic density by its c-number expecta- 
tion value = (aijaj). We will further assume that the 
bosonic atoms are heavier than the fermionic atoms as for 
a ^Li-^^Rb mixture. Since the tunneling term depends 
exponentially on the atomic mass and since the dynamic 
polarizability of the larger boson atoms is greater than 
the polarizability of the fermions leading to effectively 
deeper traps for the bosons (see Fig. 1), we can neglect 
the bosonic hopping term and retain only fermion tunnel- 
ing (23. In this limit the system can be described by an 
adiabatic Holstein model, where its ground state is found 
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adiabatically by solving the 'fast' fcrmionic problem 



(jk) j 



(2) 

treating the bosonic densities as fixed parameters and 
then adding the resulting fermion energy (which depends 
parametrically on as an effective potential to the 

'slow' bosonic Hamiltonian, 
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For Wo = 0, the adiabatic Holstein model is known to 
exhibit a Peierls instability with respect to bosonic 
collective excitations with wave number equal to Nc/M , 
corresponding to twice the Fermi wave length kp = 
Nc/2M l^^. It is favorable for the system to reduce the 
ID translation symmetry by enlarging the effective unit 
cell, thereby opening a gap in the fcrmionic spectrum. 
When the wavelength of the excitation is 2fci? this gap co- 
incides with the discontinuity of the Fermi distribution, 
so that all fermions are on the side of the spectrum which 
lowers in energy. For example, for Nc/M =1/2 the unit 
cell doubles, opening a gap in the fermionic spectrum at 
the zone boundary of the folded Brillouin zone. It should 
be noted that in difference to the standard Su-Schrieffer- 
Heeger (SSH) model [i^ used for quasi one-dimensional 
systems exhibiting a Peierls instability, the coupling to 
the bosonic degrees of freedom in our system is on-site 
whereas in the SSH model it effects the hopping proba- 
bility. 

In order to demonstrate the Peierls instability we will 
study how the energy of the system is affected by spatial 
bosonic modulations of the form 



(4) 



The density — — {jUa'-offi'^ /i)]/ gaa with /i denot- 
ing the chemical potential of the bosons, is the Thomas- 
Fermi density profile which minimizes the fixed Na 
bosonic energy H'^^j + I^Cl2j'"''j ~ ^a), in the absence 
of fermion-boson interactions. The density modulation 
depth (5n° <^ n° is generally a function of j, varying 
slowly compared to the modulation wavelength. 

Under the ansatz Q the fermion Hamiltonian (0) takes 
the form 
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where k 
and Ao 



m(a;o^Af)2/2 with m = nic ~ [gad gaa)ma, 
gac5n°j- Evidently, the mutual trapping of 
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FIG. 2: Fermionic energy spectrum for various values of the 
bosonic modulation wavenumber with luq = (a) and as a 
function of luq for a fixed modulation with k/AI ~ 1/2 (b). 



both fermions and bosons can only take place when 
gac/gaa < iTLc/nia Or fcrmion atoms will scatter out of 
the trap by the bose mean-field. In what follows we 
shall assume that this condition is satisfied. In Fig. 2 
we plot the fermionic energy spectrum, £g, as a func- 
tion of the Fermionic wave number, q, obtained from 
direct diagonalization of the fermionic Hamiltonian 
for constant Aj = A. In Fig. 2a the harmonic trap 
frequency is set equal to zero, whereas the effect of the 
trap is demonstrated in Fig. 2b by fixing the modula- 
tion frequency to k = Ad/ 2 and plotting the spectrum 
for various values of ojq. The bosonic modulation dis- 
torts the periodicity of the lattice, thereby opening a 
gap at q = fc/2. For k/M = 2kF the gap coincides 
with the Fermi momentum, so that all the states with 
\q\/M < k/2M = kp whose energy is lowered are full 
and all the states with \q\/M > k/2M = kp which in- 
crease in energy are empty. Consequently, the fermionic 
energy is minimized for k = 2kp, as depicted in Fig. El 
where we plot the fermionic ground-state energy Ec ob- 
tained by integration over the fermion spectrum up to 
the Fermi energy, as a function of the wavenumber of the 
spatial modulation in the boson field. Sharp minima are 
attained as expected, for k/M = Nc/M = 2kp. Further 
local minima of the energy, corresponding to smaller gaps 
opening at the Fermi momentum, also appear for integer 
multiples of k/M . 

The total energy of a half-filled system with k = M/2 
is plotted in Fig. as a function of the modulation 
depth A. As can be seen from the proceeding anal- 
ysis, the reduction in Ec due to the bosonic modula- 
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FIG. 3: Ground-state fermionic energy as a function of the 
modulation wavenumber for various fermion filling factors: 
Nc/M = 1/2 (solid), Nc/M = 1/4 (dashed), Nc/M = 1/8 
(dash-dotted). The external trap frequency is set to k = O.ltc 
and the bosonic amplitude modulation is set equal to tc- Ar- 
rows indicate the bosonic modulation wavenumber minimiz- 
ing Ec. 



tion is approximately linear in A. By substitution of 
Q into the bosonic Hamiltonian Q we see that the 
boson contribution to the total energy is Ea = Etf + 
{9aa/2gac) Ej ^] cos' {^Tikj / M) where Etf = (5/7)//7Va 
is the Thomas-Fermi energy. Thus, boson-boson repul- 
sive interactions increase the energy quadratically with 
the modulation depth. Consequently there is always a 
minimum in the total energy of the system Etot = Ec+Ea 
at some finite modulation amplitude indicating that the 
formation of a CDW is indeed energetically favorable. 
The optimal modulation depth decreases as tc increases 
since linear fermionic dispersion is attained at decreas- 
ingly small values of the gap. The resulting CDW can 
be detected by means of atom interferometry where the 
doubling of the unit cell will be manifested in the inverse 
lattice spacing or via the measurement of the collective 
excitation spectrum by means of Bragg spectroscopy |24| . 

Further insight into these numerical results is gained 
by employing the commonly used continuum model. For 
simplicity, we will focus in what follows, on the half fill- 
ing case Nc/M = 1/2 where the bosonic order parameter 
minimizing Ec is of the form nj=n'j+ Sn'j cos(7rj). We 
note that similar treatment can be applied for other com- 
mensurate fermion filling factors. In the continuum limit, 
the fermionic Hamiltonian (jS)) is rewritten as 

f 1 d'^ 

Hc^ dx^\x)[- — ^(J:i + A{x)a2+V{x)aQ]^{x) , 

(6) 

where 



is the spinor representation of the fermionic field in 
terms of right- and left-moving atoms, Ui are Pauli ma- 
trices ctq is the identity matrix, and m is the effective 
atomic mass. The continuum limit for the trap po- 
tential is V{x) = miJ^x^ 1^ and the gap parameter is 
Aj ^(2;)- We note that in the Takayama-Lin-Liu- 
Maki (TLM) model which is the continuum limit 
of the SSH model, there is no confining potential and 
the dispersion is linearized. Moreover (j\ appears for the 
coupling between left and right movers because the TLM 
coupling to phonons modifies the off-diagonal hopping 
rate, whereas in our case the coupling to the bose field 
modifies the diagonal self-energy terms. 

The fermion spectrum is obtained by Bogoliubov-de 
Gennes (BdG) diagonalization of the fermionic Hamilto- 
nian. We follow a similar method to the one used by 
Anderson to calculate the excitation spectrum of a su- 
perconductor with local disorder |2^. In this technique, 
which is essentially equivalent to a local density approx- 
imation (LDA), the fermionic spectrum is calculated by 
spatial averaging over spectra with different local order 
parameters. We expand the field operators "01,2(2:) as 



01 (x) 



^M„(a;)d„ , i/'2(a;) = ^w„(a;)(l„, (7) 



where dn are fermionic mode annihilation operators. The 
fimctions w„(.t), ti„(a;) are assumed to take the form 



Unix) = $„(a;)u„ , Vn(x) = $„(a;)w„, 



(8) 



where $„(a;) are harmonic oscillator eigenfunctions. Sub- 
stituting Eq. O and Eq.® into Eq. ^ and requir- 
ing that He = 'En^ndl^dn, we obtain two coupled BdG 
equations. In the framework of the LDA, we diagonalize 
the BdG equations for a local order parameter A. The 
functions ^n{x) diagonalize the spatial part of the BdG 
equations which simplify to 



(£„ - E„)Un = -lAVn 
(£„ + En)Vn = tAu„ . 



(9) 



Diagonalization of Q results in the local fermionic spec- 
trum £„ = -^/^E^^kA^ in terms of the oscillator's energy 
En ~ (jip — n + l/2)ujQ measured with respect to the 
Fermi energy, Enp = (rip + 1/2)ujq. This spectrum com- 
pares well with the numerical spectra of Fig 2b at the con- 
tinuum limit. In the limit loq 0, we have nujQ vpQ, 
where vp = 2tci is the Fermi velocity (with H set equal 
to 1) and q is the wave number for the plain wave solu- 
tion of the non-confined problem, so that the well known 
spectrum for a cos(7rj') modulation Sq = ±-\/ (qvp)^ + A^ 
is reproduced. 

Having found the fermionic local spectrum, the total 
energy functional of the system is given as the sum Etot ~ 
Ec + Ea- The fermion energy Ec is given within the LDA 
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as 



n=l n=l 

(10) 

where Aq is a constant order parameter whose value is 
the spatial average of A(a;). The boson contribution Ea 
is given by 



(11) 



where A = gac^ /{'^'^gaatc), is the dimensionless fermion- 
boson coupling constant. For a constant A(x) = Aq we 
have Ea = {gaa/'29lc)MA'^. Minimizing Etot by setting 
its variation with repect to A{x) to zero, we obtain a 
self-consistent gap equation for A(x), 
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which is similar in form to the gap equation obtained by 
Anderson |26| . 

For sufficiently wide traps, |$„(x)p can be replaced by 
its average value, thus restoring the familiar gap equation 
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FIG. 4: Total energy of the system as a function of the gap A 
for various values of tc (a) and optimal gap values as a func- 
tion of tc (b). The dashed line depicts the predicted strong 
coupling behavior of Eq. 1141 . 



1 = «fA^ dql^^ {vpq^ + ^Aj^j , (13) 

where A = 7r/2£ is a momentum cut-off of the order of 
the fermionic band width. In the weak coupling regime 
vpJ^ < Ao, Eq. lO is solved by Aq = wfAA = glJ2gaa, 
whereas in the strong coupling regime Wi?A ^ Aq we have 
the well known solution 

Ao = 2wFAexp(-l/A) . (14) 

Our numerical results agree well with these continuum 
predictions as demonstrated in Fig. 4. The weak cou- 
pling limit is confirmed by the low tc curves in Fig. 4a, 
which attain a minimum at Ao = ttA^c = I'.pAA. The 
strong coupling behavior is depicted in Fig. 4b where 
the minimum energy gap Aq is shown to precisely follow 
Eq. (|14|) (dashed line) for sufficiently large tf.. 

The formation of a CDW in the Fermi-Bose mixture 
will be affected by finite temperature effects. To observe 
Peierls instability the temperature should be lower than 
the gap T ^ Ao. The fermionic excitation spectrum de- 
pends exponentially on the ratio T / Aq so that the num- 
ber of excited fermionic atoms should be exponentially 
small. However the density of thermal bosonic excita- 

1/2 

tions has a power-law behavior Uex = {Na/V){T /Tc) , 
where V is the volume and Tc is the critical tempera- 
ture. This leads to a more restrictive constraint T /Tc <C 
(Ao/gac'T-)^, where n is the bosonic density. 
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